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SYMBOLS 


Symbol Name Units 

K. E. Kinetic energy J 

Kinetic energy per unit of 
rotor weight, termed energy 

densily I/N 

a Maximum centrifugal stress Pa 

Radial centrifugal stress Pa 

Circumferential centrifugal stress Pa 

Kj Flywheel rotor shape factor — 

E Kinetic energy per unit of swept 

vs ^ n 

rotor volume 

W/V^ Rotor weight per unit of rotor N/m^ 

swept volume 

T Thickness of rotor or of a square single 

filament 

R Radial thickness of a round single 

filament 

7 Weight density of material N/m^ 

a Inner rotor radius ni 

b Outer rotor radius ni 

V Poisson ratio of contraction in the radial 

direction due to extension in the circumfer- 
ential direction 

ID Inside diameter m 

OD Outside diameter ni 

CO Angular velocity rad/s 

N Rotations per minute rot/min 
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Symbol 


Name 


Units 


W(co) 

Parasitic losses in an energy storage system— 
a function of angular rotor speed 

W 

P(co) 

Load or input power for an energy storage system— 
a function of angular rotor speed 

W 

r(w) 

Efficiency of power transfer in an energy 
storage system— a function of angular rotor speed 

- 

jg 

Flywheel polar moment of inertia 

Js^ 


Total load on an energy storage 
system W (co) + P (co) 

W 

n 

The number of single filaments in a “brush” 
superflywheel 

- 

P (laminar) 

Windage losses assuming laminar fluid 
flow 

W 

P (molecular 
flow) 

Windage losses assuming molecular stream 
fluid flow 

w 

r 

Radius from center of rotation 

m 

Wh 

Watt hours 

Wh 

Ibf 

Pounds force 

Ibf 
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INTRODUCTION 

One of the main purposes of energy storage is to conserve petroleum resources through the 
reduction in both the dollar and the environmental costs of providing energy to the con- 
sumer. The basic problem facing utility companies is, and will continue to be, variations 
in the demand for electric power. Kalhammer and Zygielbaum (reference 1), Keller 
(reference 2), and Rabenhorst* point out that electric utilities use: 

• Base load generators to serve the load which continues 24 hours a day. This 
function is served by the more efficient fossil fuel or nuclear generators which 
operate nearly 100 percent of the time at full throttle. Typically, the base 
load is approximately 45 percent of the peak load and 70 percent of the 
system power output. 

• Intermediate load generators serve the broad daily demand peak. This 
function is served with older, less efficient fossil fuel steam generators and 
gas turbines; these are normally shut down at night. 

• Peak load generators serve the immediate load peaks of the day. These are 
typically gas turbines or diesel engines. 

As Keller points out, the base load generators may supply about 70 percent of the total 
systems power at the lowest delivered cost, while the intermediate load generators will 
furnish about 25 percent of the power at significantly higher cost. The remaining power 
is derived from peak load generators which, although relatively inexpensive in terms of 
capital investment, require costly and scarce special fuels (such as number 2 fuel oil), 
operate at low efficiencies, and necessitate high expenditures for maintenance. 


*Rabenhoist, D.W., “Use of Flywheels for Energy Storage,” presented at American Chemical Society, Fall 1974 
Meeting of Energy Storage Symposium, Atlantic City, New Jersey, September 12, 1974. 



With the cost of fuel increasing rapidly, it is extremely expensive to run the intermediate 
and peak load generators. Laaspere and Converse (reference 3), for example, point out that 
the Central Vermont Public Service Corporation has introduced a rate option where electri- 
city used between 8 and 1 1 a.m. and between 5 and 9 p.m. costs six times more than that 
used at other times. This type of costing differential provides the economic impetus for 
the development and installation of high efficiency energy storage systems. 

In the past, electric utilities have attempted to reduce the high cost of peak power genera- 
tion through pumped hydroelectric storage (reference 4). Typically, water is pumped up 
hill during off peak hours and stored in a reservoir for later peak hour use. This method is 
approximately 66 percent efficient (reference 4) and is limited to those locations where 
geography permits the construction of large reservoirs. As an example, a 1900-MW, 
15,000-MW h system is jointly operated by the Consumers Power Company and Detroit 
Edison Company. The storage reservoir is a manmade lake more than 3 by 1.5 km (1.86 
by 0.93 mi) constructed at a cost of more than $300 million. To deliver the required 
amount of power, roughly 10® m® of water must be pumped and discharged, and the overall 
system represents a capital cost of approximately $200/kW. Furthermore, the stored energy 
density is quite low (reference 2)— 81 J/N (0.1 W h/lbf). 

There are numerous mechanisms which are theoretically capable of storing energy. The 
following list describes some of those mechanisms which are currently under consideration 
for electric utility storage systems; 

• Mechanical— flywheels (kinetic energy), 

• Mechanical— pumped water storage (potential energy), 

• Mechanical— compressed air storage (potential energy), 

• Chemical— fuel cells, 

• Electrical— secondary storage batteries, and 

• Electrical— superconducting magnet. 

Not all of these mechanisms are sufficiently technically developed nor economically feasible 
for application to a complete system of electric utility energy storage. 

Kalhammer and Zygielbaum (reference 1) have compared various storage systems which 
utilize the six mechanisms previously mentioned, and they conclude that flywheel energy 
storage will be viable in utility application when acceptable reliability and long lifetimes 
can be assured. They suggest that lifetimes greater than 20 years and power conversion 
efficiencies greater than 60 percent will be needed to make flywheels applicable to peak 
shaving for distributed utility energy storage systems. 

The concept of energy storage in rotating flywheels is not new; Clerk (reference 5) has 
provided a useful history of the development of flywheel systems. Clerk points out that 
the most severe practical limitations of stored kinetic energy systems have been parasitic 
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losses due to friction, windage, relatively poor energy storage per unit weight of flywheel, 
and the hazard of catastrophic failure of the rotor. 

Before discussing recent work on flywheel energy storage systems it is useful to consider 
measures of performance which are used to compare maximum storage capabilities among 
various flywheel configurations. The following quantities have been used by numerous authors 
throughout the literature. 

Symbol 

Description 

Unit 

E 

W 

Energy density— stored kinetic energy per 
unit of flywheel weight 

J/N (W h/lbf) 

Evs 

Stored kinetic energy per unit of swept fly- 
wheel volume 

Pa (psi) 

W/V 

Flywheel weight per unit of flywheel swept 
volume 

N/m^ (Ibf/in.^) 

Ks 

Nondimensional factor for comparing the 
effectiveness of various flywheel shapes— 
maximum 1.0 

- 

oh 

Flywheel specific strength, maximum working 
stress in flywheel per unit of weight density 

J/N (W h/lbO 


The following equation (discussed in the section on flywheel rotors) governs the energy 
density for any given flywheel configuration. 


Also, the following relationship governs and W/Vs. 

W/Vs = E^^/E^ (2) 


Useful conversion factors for manipulating many of the energy storage units found in the 
literature to a common system are shown in appendix A. 

Rabenhorst (references 6 through 13) was among the first investigators to show that 
anisotropic (composite) materials have capabilities of extremely large specific strengths in 
flywheel applications. For example (reference 6), a 70 percent graphite whisker/epoxy 
has a usable specific strength of 457 kJ/N (565 W h/lbf), while the material Kevlar-49 
(reference 13) has capabilities of 88 kJ/N (109 W h/lbf). These values represent a 
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substantial improvement over maraging steel at 26 kJ/N (32 W h/lbf) and at least an order 
of magnitude over the lead -acid storage battery (reference 6) at 6.4 kJ/N (8 W h/lbf). 

Because of the large values of specific strength for anisotropic materials, those kinetic storage 
systems which allow optimal use of essentially straight anisotropic materials are termed 
superfly wheels. Rotor geometry for composite materials, however, is limited to those geo- 
metries that develop nearly unidirectional centrifugal stresses. 

There are two basic rotational configurations which develop essentially uniaxial stress pat- 
terns, allowing maximum use of the composite material specific strength. These are; 

• The fanned circular brush superflywheel as described by Rabenhorst (references 
8 through 10). In this configuration a number of thin anisotropic rods are 
spinning about their minor axis (K^ = 0.33 for this configuration). 

• The thin rim or multirim super fly wheel, as described by Rabenhorst (reference 
11) and Post et al. (reference 14). In this configuration a number of independent 
thin rings (spinning about a common axis) are essentially in pure circumferential 
uniaxial stress (K^ = 0.50 for this configuration). 

If either of these flywheel configurations is used in an energy storage system, there must be 
rotation about a common shaft and either a connection hub (for the brush) or spokes (for 
the multirim) are required to connect the rotor to the power shaft. Invariably, the super- 
flywheels must operate in a vacuum environment, forcing the designer of a total storage 
system to deal with potential problems in such items as seals, bearings, pumpdown equip- 
ment, enclosures, and electrical/mechanical converters. The magnetically levitated spoke- 
less ring flywheel system described in this report overcomes a number of these potential 
problems. 

As of this writing no complete anisotropic flywheel storage system has been built and tested, 
although Dugger et al. (references 15 and 16) have evaluated individual single filaments of 
composite materials in spin tests. Interestingly however, a small 222-N (50-lbf) magnetically 
suspended, anisotropic ring for angular momentum spacecraft control is currently under 
test at Langley Research Center (reference 17). A number of the advanced concepts used in 
this ring were developed at the Goddard Space Flight Center (references 18 through 20) and 
the application of these concepts to the levitated ring energy storage systems will be described 
later. 

The purpose of the remainder of this report is to; 

• Develop a working equation and set general requirements for the kinetic 
energy storage system, 

• Compare rotor geometries for various flywheel configurations. 
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o Discuss the magnetically levitated spokeless ring in terms of the required 

modules which are necessary if the concept is to be implemented as an energy 
storage system, and 

• To suggest a design for a 7 kW h magnetically levitated spokeless ring system 
which can be used to demonstrate the system. 

THE GENERAL ENERGY STORAGE SYSTEM 

A block diagram of a rotational kinetic energy storage system is shown in figure 1 . The 
losses of both the storage system and of the electrical conversion elements, which are con- 
stantly present parasitic losses, are shown as watts. In general, these losses are a function 
of the flywheel angular rotational frequency (co) and are present during both charge and 
discharge cycling. 



Figure 1. Block diagram of flywheel energy storage system. 


As soon as power flow begins (either charge or discharge), the efficiency of the converting 
elements (which can also be a function of angular rotational frequency) enters the calcula- 
tion. 

If an equation for the energy balance of this system is written, then: 


-W(w)-P(co)/f(w) 

Discharge cycle 

(3) 

-W(w) + P(cu)/r(oj) 

Charge cycle 

(3a) 


Units 


where: 

K.E. 

= stored kinetic energy 

joules 


t 

= time 

seconds 


W (co) 

= losses, a function of co 

watts 
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+ P (co) = power input (charge), a function of co watts 

- P (co) = load output (discharge), a function of co watts 

f (co) = efficiency of power conversion, a function 

of CO 

CO = angular rotational frequency radians/second 

If the moment of inertia of the flywheel is constant, then (3) may be rewritten as: 

^ ^ -WM P(^) ( 4 ) 

dt Jg Jg? (co) 

Units 


where: co = angular rotational frequency radians/second 


Jg 

= flywheel polar moment of inertia 

joules • second 

For the simple case where the load and losses are independent of co, the kinetic energy de- 
creases linearly with time while co decreases as the square root of time. For this case: 


9.14X10-^Jg 

t* = (N2 - 

(5) 



Units 

where: t* 

= time to slow to N^. 

minutes 

N. 

= initial rotational speed 

revolutions 
per minute 

Nc 

= final rotational speed 

revolutions 
per minute 

Jg 

= flywheel polar mass moment of inertia 

joules • second^ 

Wt 

= total load 

watts 

Those variations of W, P, or f (as functions of co) which are not easily expressible mathe- 
matically would require numerical solution of equation (4). This procedure would be used 
to generate co versus t plots which could then be converted directly to K.E. versus t (since 
K.E. = V 2 Jg co^ ) plots and the time for the system to drop to some percentage of the initial 
kinetic energy could be obtained directly from the plot. 
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Based in part upon equations (1) and (5), the requirements of a flywheel storage system are: 

• Efficient rotor shapes to take full advantage of rotor strengths, 

• Lx)w parasitic losses, 

• Reliable operation over long lifetimes, 

• High efficiencies of power conversion, 

• Minimization of damage in the event of rotor failure, and 

• Ability to avoid self-destruction in the event of seismic disturbances. 

FLYWHEEL ROTORS 

The study of flywheel rotor shapes has received the greatest attention of all the components 
of the flywheel energy storage system. Lawson (reference 21) and Dugger et al. (references 
15 and 16) have analyzed the stored kinetic energy in various flywheel shapes (where cen- 
trifugal stresses limit performance) and have produced a listing of numerical shape factors 
(K^ in equation (1 )) for particular geometries. Unfortunately, the numerical shape factors 
depend on geometrical properties of the specific rotor shape, and general formulas are not 
provided by either Lawson or Dugger. 

The problem of rotor shape characterization is best visualized by referring to figure 2. For 
any rotor shape, the following information should be expressible in mathematical form in 
order to evaluate rotor performance: 

• K.E. — stored kinetic energy in the rotor, 

• E^ — kinetic energy per unit of rotor weight, 

• ^max ” maximum (limiting) allowable centrifugal stress in the rotor, 

• Kg — the nondimensional shape factor for a given geometry, 

• E^g — kinetic energy per unit of swept volume, and 

• W/Vg — rotor weight per unit of swept volume. 

For all of the assumed flywheel geometries which are analyzed in this report, the following 
procedure is used to characterize the rotor: 

1 . An expression is written for the kinetic energy of the rotor, 

2. The kinetic energy expression is divided by the weight of the rotor to 
form E . 

w 
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w ^ ^ 


ROTOR 


Figure 2. Flywheel rotor characterization. 

3. An expression is written for the magnitude of the maximum centrifugal stress, 
and this is used to eliminate cj from the expression. 

4. The expression for is multiplied by the rotor weight and divided by the 
swept volume to obtain an expression for E^^. 

5. The expression for W/V^ is obtained by dividing E^ into E^^. 

Applying steps 1 through 3 results in an expression of the form: 

E = K X -2_ 

W s ^ 

which has been given previously in equation (1). In this expression, depends only on 
rotor geometry. 

General rotor shape characterization formulas are shown in appendix B for the following 
geometries: 

1. Constant stress disk— A rotor shape selected such that a = a. = = constant, 

r C/ iTlaX 

everywhere in the rotor (reference 22). 

2. Pierced disk (ring)— A constant thickness rotor analyzed for both radial and 
circumferential stresses (treated as isotropic for this analysis). 

3. Hoop— A thin (ID/OD approaches 1.0) constant thickness pierced disk, where it 
is assumed = 0. A given swept volume could be occupied by a number of in- 
dependent hoops with the effect of the overall ID/OD ratio being much less than 
1 .0. For this case the hoop formulas would still apply. (This was the multirim 
analysis used by Rabenhorst (reference 1 1).) 

4. Round single filament— A single round rod rotating about its minor axis. 

5. Square single filament— A single square rod rotating about its minor axis. 

6. Brush— n single filaments (either round or square) all rotating about a common 
minor axis and in the same thickness plane (references 8 and 13). 
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The particular geometries were selected because of the following reasons: 

• They are the only configurations which provide nearly uniaxial centrifugal 
loading, thus allowing the use of the maximum strength properties of anisot- 
ropic materials. 

$ They have been suggested by both Rabenhorst and Post as “best” rotor geometries. 

A comparison of the general formulas for the above six shapes is shown in table 1. 

To compare the shape factors of various rotor geometries, specific inside-to-outside radius 
a/b must be selected. Based upon Rabenhorst’s discussion of the brush and multirim 
flywheels (reference 11), two cases are suggested. 

1. The a/b ratio where the volume occupied by the rotor divided by the volume 
swept by the rotor is the same as for the brush flywheel (approximately 15 
percent). This requires a/b = 0,92. 

2. The a/b ratio where the energy density for the hoop is the same as the energy 
density for the single filament and brush. This requires a/b = 0.58. 

A comparison of rotor performance for these two cases is shown in table 2. For the square 
single filament and brush, it has been assumed n = 20, r = 1 .2 cm (0.5 in.) and 2b = 2.25 m 
(88 in.) after Rabenhorst and Taylor (reference 13). 

Table 1 

Comparison of Rotors 


Shape 


E^ 

w/v, 

Constant sties'^ Jisk 
ID = 0 

OD = oo (best case) 

1.000 X - 

X a 

7 

Pic iced disk 

(h) va 


H’) 

A ^ 

3 + y + ( ^ y ( 1 - y) 

Hoop 

lb/ y 2 

4 7 


X 
1 

Single filament 
Square 

Round 


K X o X 

7Tb 



K X.X (K) 

"{ 1 ) 

Brush square single 
filaments 

1x2 
3 7 

K X 

* 7T b 

V 

7X 

7T b 
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For electric utility storage systems, long life, reliable operation, high efficiency, and low 
cost are primary requirements. However, only the rotor is being compared in table 2. Fur- 
thermore, the comparison does not include fabrication costs or other strength limiting 
factors (besides centrifugal stresses) which may arise when a rotor is assembled as part of a 
complete storage system. Bearing in mind these limitations, the following conclusions are 
drawn: 

• Based on , the hoop and pierced disk are better than the brush for inside- 
to-outside diameter ratios > 0.58. 

• The hoop and pierced disk are nearly identical for both , E^^, and W/V^. This 
is expected because the hoop was analyzed as a special case of the pierced disk. 

• The constant stress disk is limited to isotropic materials because = Og = a = 
constant throughout the disk. 

• For the pierced disk or hoop to be competitive with a constant stress disk, the 
composite material must have a usable specific strength (a/7) greater than 2.13 
(1.00/0.47) of the best isotropic material which is potentially available for a 
constant stress disk. Actually, this is a worst case condition since, as Lawson 
(reference 21) points out, the best practical constant stress disk is limited to 

=0.93. 

A material comparison listing taken from various references as indicated is shown in table 3. 
Assuming that maraging steel is the highest strength isotropic material for a constant stress 
rotor (reference 21), then several composite materials have improvement factors greater 
than 2. In particular, Kevlar-49 is a prime candidate for a composite rotor (it is now selling 
for approximately $2 to $2.50/N ($8 to $10 /Ibf) but the future price may drop to about 
$0.50 to $0.70/N ($2 to $3/lbf) (reference 12)). 

A study of energy storage rotors for power and attitude control of spacecraft, an Integrated 
Power/ Attitude Control System (IPACS) has been reported by Notti et al. (references 23 and 
24) and Keckler and Jacobs (reference 25). The authors considered the use of the three 
types of composite materials: 

1. Kevlar-49 (PDR -49-III) 

^ max = 1 .38 G Pa (200,000 psi) 

7 = 13.6 kN/m^ (0.05 Ibf/in.^) 

2 = 102kJ/N(126Wh/lbO 
7 
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Table 2 

Comparison of Rotor Performance 


Shape 


E, 

w/v^ 

Constant stress 
disk (O.D. «>) 

1,000^ 

7 

1.00 a 

1.00 7 

Pierced disk (ring) 
(7 = 0.3) 

r = 0.92: 0.47 

b 7 

0.07 a* 

0.157 


^ = 0.58:0.38 ^ + 

b 7 

0.25 a+ 

0.66 7 


^ = 0.001 ;0.30S 
b 7 

0.30 a 

1.00 7 

Hoop 

{ 1 or many) 

= 0.92; 0.46 

b 7 

0.07 0 * 

0.15 7 


^ =0.58; 0.33 
b 7 

0.22 a+ 

0.667 

Single filament 
square 

0332: 

7 

0.002 a 

0.007 7 

Brush 20 siiLiare single 
filaments 

0.33 S 
7 

0.05 a 

0.15 7 


“ 0 92 I Ills tor responds lo W/V (or ihe pierced disk and hoop = W/V for the brush. 


+ - - 0.58 

b 


Tins corresponds to for the hoop = for single filament and bru.sh. 


2. Graphite Epoxy 

° = 1.24 G Pa (180,000 psi) 

7 =15.7 kN/m^ (0.058 Ibf/in.^) 

a =78.8 kJ/N (97.4 Wh/lbf) 

7 

3. Boron Fibers 

a = 2.41 G Pa (350,000 psi) 

7 =25.5 kN/m^ (0.094 Ibf/in.^ ) 

a = 94.7 kJ/N (117 Wh/lbf) 

7 

in a variety of geometrical configurations. 

A listing of the results of their study is shown in table 4-A, This may be compared to the 
results shown in table 4-B which were prepared for a ring with ID/OD = 0.9 (K^ = 0.45). Note 
that the specific energy of the ring rotors are greater than all entries in table 4-A. 
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Table 3 

Comparison of Material Properties 


Material 

Ultimate Tensil Strength 
(UTS) 

N/m^ X 10® (Ibf/in.^ X 10®) 

Weight Density D 
kg/m® X 10® (Ibf/in.®) 

Usable Strength 
S = UTS X 72% 

N/m^ X 10® (Ibf/in.^ X 10^) 

S/D 

J/N X 10® (Wh/!bf) 

Improvement 
Over Maraging 
Steel 

Music Wire* 

41.34 (600) 

78.10 (0.283) 

43.48 (431) 

0.39 (48) 

1.5 

Boron/M agnesium * 

10.82 (157) 

24.01 (0.087) 

9.30 (135) 

0.39 (48) 

1.5 

Steel/Epoxy* 

29.07 (422) 

48.71 (0.1765) 

21.22 (308) 

1 0.44 (54) 

1.7 

Future Music Wire* 

68.90 (1000) 

78.11 (0.283) 

49.61 (720) 

0.65 (80) 

2.5 

Beryllium Wire* 

17.23 (250) 

18.49 (0.067) 

12.40 (180) 

1 0.68 (84) 

2.6 

Future Boron/Magnesium* 

24.12 (350) 

24.29 (0.088) 

17.36 (252) 

0.73 (90) 

2.8 

Boron/Epoxy* 

22.03 (320) 

' 20.42 (0.074) 

15.85 (230) 

' 0.79 (98) 

3.0 

HTS Glass/Epoxy* 

22.74 (330) 

20.53 (0.0744) 

16.39 (238) 

0.75 (93) 

2.9 

1 Present Graphite/Epoxy* 

19.50 (283) 

16.00 (0.058) 

14.12 (205) 

0.88(109) 

3.4 

j Future Glass/ Epoxy Prestressed* 

32.38 (470) 

20.42 (0.074) 

23.29 (338) 

1.16(143) 

4.5 

Future Graphite Epoxy* 

34.45 (500) 

16.56 (0.060) 

24.80 (360) 

1.52(188) 

5.8 

Theoretical Future* 

84.75 (1230) 

17.94 (0.065) 

60.98 (885) 

3.73(461) 

14.3 

70% Graphite Whisker/ Epoxy* 

103.35 (1500) 

16.56 (0.060) 

74.41 (1080) 

4.57(565) 

17.6 

Kevlar-49t 

16.54 (240) 

13.80 (0.05) 

19.92 (173) 

0.88(109) 

3.4 

Maraging Steel§ 

28.18 (409) 

79.76 (0.289) 

20.26 (294) 

0.26 (32) 

1.0 


*See reference 6. 
tSee reference 13, 
§See reference 21. 



Table 4 


Comparison of Rotor Shapes and Performance 


A. Comparison of Rotor Shapes (Reference 25) 

Design 

Material 

Specific Energy 
W h/kg (W h/lb) 

Circular brush 

Boron 

67.75 (30.76) 

Orthrotropic disk 

Graphite 

28.24 (12.82) 

Tape wound 

PDR-49-III 

71.10 (32.28) 


Graphite 

60.37 (27.41) 

Isotropic 

Titanium 6A1-4V 

47.75 (21.68) 

Isotropic 

Steel 300M 

34.34 (15.59) 

B. Ring Rotor Performance with ID/OD ^ 

= 0.9 (K^ = 0.45) 

Ring 

Boron 

116.08 (52.7) 

Ring 

Graphite 

96.47 (43.8) 

Ring 

PDP-49-III 

124.89 (56.7) 


One entry of interest in table 4-A is the orthrotropic disk (that is, pierced disk) which has 
an effective shape factor of only 0.13 (12.82/97.4). The reason this is much lower than the 
corresponding entry in table 4-B (for the graphite ring) is that the shape of the orthotropic 
disk is not thin, and the ultimate speed is limited by the radial delaminating centrifugal 
stress— a problem not encountered for the ID/OD = 0.9 ring. However, the authors feel that 
a major reason a thin ring was not considered in the IP ACS study was because there was no 
easy way to couple ring motion to a power source without derating the specific energy of 
the ring. A method to overcome this difficulty is suggested later in this report. 

The results of the rotor study may be summarized as follows: 

• The stored energy per unit of weight makes a composite rotor ring (that is, 
pierced disk with a/b = 0.9) flywheel attractive as an energy storage device. 

Its energy density is superior to a brush rotor for all practical ID/OD geometric 
ratios. 

• The principles of magnetic suspension and brushless commutation (references 
26 and 27) can be applied to energy storage rings just as they have been to 
angular momentum control rings in spacecraft (references 17 through 20). 

• The feasibility of constructing a composite storage ring has already been 
successfully demonstrated for a 222-N (50-lbO rotor (reference 17). 


13 



THE RING FLYWHEEL STORAGE SYSTEM 


Figure 3 is a block diagram of the major modules of the complete ring storage system. It is 
suggested that the storage ring be magnetically suspended and enclosed in a toroidal vacuum 
container with all the suspension and drive electronics completely enclosed— only wires 
pierce the enclosure. 

The storage module would consist of a flywheel ring rotor configuration as shown in figure 
4. The ring could be constructed of a high strength composite material (for example, 
Kevlar-49) and contain the magnetic suspension components on the inner radius. The 
center of the inner radius will be used for the motor/generator salient poles, which will be 
attached to the rotor. The ring will carry only iron material, and only the inner radius of 
the ring will be used for this purpose. 


SUSPENSION ENVIRONMENT 



~I 

I 


L 
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Figure 3. Block diagram of the ring flywheel energy storage modules. 

Because the ring is now composed of both a composite material and iron, the overall shape 
factor will be less than for a single composite material ring. However, the reduction in shape 
factor may be minimized by proper selection of the magnetic iron material (i.e., trading off 
some “magnetic softness” for an increase in mechanical strength) or, through an interfer- 
ence fit between the iron and composite ring or, through the use of new high mechanical 
strength magnetically soft “Metglas” (reference 28). 

The environment module would consist of a toroidal enclosure which would house both the 
rotating ring and the levitation and motor/generator components. This module would con- 
sist of vacuum pumpdown equipment which will allow the enclosure to be placed into the 
0.13 to 0.0013 Pa (10"^ torr to 10’^ torr) pressure range. At these pressure levels the parasitic 
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windage losses will be negligible, and this configuration will not require seals. It is also en- 
visioned that the continuous pumpdown requirements will be minimal since the integrity 
of the enclosure can be easily maintained. 

The suspension module will consist of a circumferential stationary force coil and an iron 
ring attached to the rotor. The magnetic circuit will establish a magnetic flux across the 
gap between the moving ring and the stationary coils. Control of rotor position will be 
passive in the axial direction and active in the radial direction. Because of a unique per- 
manent magnet flux biasing system, the power requirements for radial control will be near 
zero.* 

The motor/generator of the input/output module will consist of a noncontacting, electrically 
commutated stationary armature which cuts the flux lines of a salient pole variable 
strength moving field. The electrical output will consist of a periodic voltage (constant 
amplitude-variable frequency) which will be rectified and then inverted to 60 Hz, 1 1 0 V 


*Simpson, P.A., “Annular Momentum Control Device (AMCD) -Magnetic Suspension System,” Final Report, Sub- 
contract Number 01397, Prime Contract Number NASl-12529, Cambridge Thermionic Corporation, Cambridge, 
Massachusetts, 1974. 
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power. To make the generator function as a motor, 60 Hz 1 10 V a.c. power will be rectified 
and fed to the electronic commutation network. This will produce an accelerating torque on 
the ring which will bring it up to speed. 

A comparison of the system requirements for a brush or multirim “conventional” flywheel 
system and the spokeless levitated ring system are shown in table 5. The primary advantages 
of the levitated ring concept are: 

• Elimination of all conventional bearings (and associated lubrication problems) 
resulting in increased reliability and long life; 

• Elimination of vacuum seals (only electrical wires will penetrate the enclosure 
system)— also, the need for large continuous pumpdown systems (due to seal 
leakage) is reduced if not completely eliminated (if outgassing is eliminated a 
one shot pumpdown will suffice); and 

• Elimination of all mechanical contact by the input/output electronically commu- 
tated armature resulting in improved reliability and long life. The electrical 
output will be in the form of alternating voltage and the output from individual, 
reasonably sized, independent systems can then be combined through a trans- 
former. 

DISCUSSION OF THE RING STORAGE MODULES 
Storage— The Ring Rotor 

The ring rotor (refer to figure 4) contains the energy storage mass (composite ring) along 
with the “magnetically soft” iron which is required for both the magnetic suspension (sus- 
pension module) and the motor/generator (input/output module). The presence of iron 
will cause an effective reduction in the shape factor of the total ring, but this effect can 
be minimized (and the overall shape factor caused to approach that of a composite material 
ring acting alone) by trading off “magnetic softness” with increasing alloy content to im- 
prove the mechanical strength of the iron. For purposes of this report, it is assumed that the 
presence of all iron is negligible (there will be no speed derating due to the iron), and the 
rotor will be analyzed as if it were simply a constant thickness composite material rotor. If 
this is true, then the ring rotor must operate under the following strength constraints: 

• The speed of the rotor must not cause the maximum circumferential stress 
to exceed the working stress of the composite material. 

• The speed of the rotor must not cause the radial (delaminating) stress to 
exceed the working cross-fiber strength of the composite material. 

Morganthaler and Bonk (reference 29) have shown that the radial and circumferential 
stress equations for an anisotropic pierced disk approach the classical isotropic formulas 
(reference 22) as the ID/OD ratio of the disk exceeds 0.9. 
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Table 5 

Comparison of System Requirements for a Conventional versus 
Spokeless Levitated Ring Flywheel System 


Conventional Flywheel System 

Levitated Ring 

ROTOR 


Hub 

Single ring storage 
element (with imbedded 

Shaft 

Storage elements (multirim or 
brush) 

Shaft/storage element connection 
spokes 

iron material) 

SUSPENSION 


Conventional ball bearings 

Magnetic levitation 
(passive axially, active 
radially) 

ENVIRONMENT 


Vacuum operation 

Vacuum operation (only 
1 pumpdown to 0.13 to 

Seals required 

Continuous pumpdown may be 
required at less than 1 .3 N/m^ 
(10“2 torr) 

0.0013 N/m^ (10“^ to lO'" 
torr) required) 

Enclosure system must be larger 

Enclosure system is a 

than OD of flywheel 

toroid with unenclosed 
center 

MOTOR/GENERATOR 


Synchronous converter 

Synchronous converter not 

required to hold output voltage 

required, special winding 

constant as shaft speed changes 

provides constant voltage 
as speed changes 

Rectifier/inverter required 

Rectifier/inverter required 
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The flywheel ring may thus be analyzed using the conventional isotropic stress formula pro- 
vided the above condition is shown to be true. Therefore the following conditions are 
assumed : 

• The classical isotropic and Og formulas apply to the composite ring rotor. 

• The ring will not delaminate as long as the maximum transverse stress 
(Oj.) is held below 0.5 percent of the maximum fiber stress (Og ). Accord- 
ing to Post et al. (reference 14), this is a very reasonable assumption. 

Thus it can be shown (appendix C) that the ring will not delaminate provided the ID/OD 
ratio is greater than 0.90). 

A greatly simplified design outline for sizing a composite ring rotor is shown in appendix D. 

It has been included to provide an estimate of rotor size for a 7-kW h system so that typical 
numbers (for parasitic losses and suspension and motor/generator shares) may be obtained 
in the sections which follow. The results for a 1 0-kW h storage ring are summarized in table 6. 

Test results by Rabenhorst and by Dugger et al. have shown that when composite fiber rods 
fail they do so by delamination instead of rapid crack propagation. Typically, only 
2 percent of the rotational kinetic energy was actually transported to the rotor container 
during overspeed failures. This may be contrasted with steel rotor failures which absorb 
little energy and cause the rotor to break into as few as three large pieces. It is thus antici- 
pated that the containment problems of an overspeed failure for a ring composite rotor would 
be minimal. Furthermore, the ring system described in this report could easily be placed 
underground, completely eliminating any hazard to human life in the event of a failure. 

Table 6 

(7 kW h full to half speed) 

Shape and Speed for a 1 0-kW h Storage Ring made of Kevlar-49 Composite Material 


Weight 

1085 N (244 Ibf) 

OD 

2.0 m (79 in.) 

ID 

1.8 m (71 in.) 

Thickness 

13.4 cm (5.3 in.) 

Angular speed 

855 rad/s (8163 RPM) 

OD tip speed 

855 m/s (2805 ft/s) 

Mach tip speed 

2.54 
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Environment 

The purpose of controlling the ring environment is to minimize windage losses of the spin- 
ning rotor. As the tip speed of the rotor increases and the enclosure pressure decreases, the 
fluid begins to look increasingly like a molecular stream, and the conventional drag co- 
efficients (based only on the Reynolds number) no longer apply. Rohsenow and Choi 
(reference 30) provide drag coefficient formulas for the molecular stream region, and these 
have been incorporated into the windage loss equations which are derived in appendix E. 

A summary of the total windage losses for the rotor described in table 6 is shown in table 7. 
The losses P (laminar) are those which are calculated by assuming that the flow remains 
laminar as pressure decreases; losses are overestimated considerably at Pa (10'^ torr). 

Based on the results shown in table 7, it is reasonable to assume that the parasitic windage 
losses are negligible at pressure of 0.0013 Pa (10"^ torr). This magnitude of pressure is 
easily obtained with conventional pumpdown equipment, and maintaining it in a com- 
pletely sealed enclosure will not be difficult. 


Table 7 

Total Windage Losses for 7-kW h Rotor Described in Table 4 


Pressure 

P (laminar) 

P (molecular flow) 

Pa (torr) 

(watts) 

(watts) 

0.13 (10-^) 

193 

225 

0.013 (10^) 

61 

31 

0.0013 (10-®) 

21 
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Suspension System 


The purpose of the suspension system is to provide stable support for the ring rotor in both 
the axial and radial directions with the minimum amount of power loss (that is, minimum 
parasitic loss) and longest possible mechanical life. 

Based upon these requirements, a modified electromagnetic system was chosen. A cross 
section schematic of one of the two suspension coils, previously shown in figure 4, is 
shown in figure 5. The principle of operation of this system is to establish a steady-state 
magnetic flux in the air gap (by the samarium-cobalt permanent magnets) and allow the 
change in magnetic reluctance, due to axial sag, to support the rotor weight without active 
axial control. The exact amount of axial sag can be reduced to any desired value by simply 
increasing the number of teeth on the moving iron. If the radial gap clearance is sensed, 
then the radial rotor position can be actively controlled by varying the current through the 
stationary force coil, thus holding the radial gap clearance constant. (At least one direction 
must be active for the system to be stable.) 
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PERMANENT MAGNET 
(continuous around inside diameter) 


Figure 5. Cross section of magnetic suspension. 


Shown in appendix F are the detailed design equations for the 1 0-kW h rotor described in 
table 6. These results are summarized in the following list. 


Rotor Suspension Characteristics 


Axial Radial 


Passive, Minimum Reluctance 

Servoed 


Max. Load 

= 222 kg (488 lb) 

Nominal Flux 
Density 

= 0.6 1 5 tesla (40kl/in^ ) 

Rotor ( ID) 

= 1.8 m (70.8 in.) 

Control Flux 
Density 

= ± 0. 1 54 tesla (± 1 0 kl/in^ ) 

Clearance 

= 0.152 cm (0.06 in.) radial 

Peak Force 

= 11 10 kg (2440 lb) 

No. of Teeth 

= 16 

Control 

Sensitivity 

= 3.36 kg/W (7.4 Ib/W) 

Tooth Width 
Stiffness 

= 0.1 57 cm (0.062 in.) 

= 6930 kg/cm (8000 Ib/in.) 

21 Segments 
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The type of electromagnetic forcer described in this report has been successfully applied to a 
223-N (50-lbf) ring rotor intended for angular momentum control in spacecraft applications* 
(reference 17). It is anticipated that the radial gap sensors and the force coil electronics 
(to provide the current for the stationary force coils) will be directly applicable to the spoke- 
less ring storage system described in this report. 

It is further anticipated that the parasitic wattage losses of the suspension system will be 
small. Eddy current losses, due to the movement of magnetic field, will be eliminated 
by making the stationary force coil continuous around the inner circumference of the ring. 
Hysteresis losses, due to the cycling of a magnetic field in the iron, will be dependent on the 
amount of mechanical tolerance in the shape of the ring and the perturbing forces which 
act on the storage system. Because it is doubtful that either of these factors will be large, 
the hysteresis losses are considered negligible. The only other loss is the quiescent power 
dissipation of the control circuit for the force coil. Based on a spacecraft momentum con- 
trol device, it is felt these losses can be made quite small (a few watts). 

An additional advantage of this magnetic suspension system is that there will be little 
penalty in hysteresis losses if it is necessary to use a stronger iron alloy for the moving iron. 
The fact that a stronger alloy can be used provides a firmer basis for the previous assumption 
of treating the storage ring as a single element of composite material (no iron considered). 

Motor/Generator 

The purpose of the motor/generator is to provide for the input and output of electrical 
energy from the moving rotor. The general requirements for the motor/ generator system 
are: 

• Be compatible with the spokeless magnetically levitated ring geometry, 

• Provide long life and high energy conversion efficiencies, 

• Have negligible parasitic losses, 

• Be functional as both a motor or generator, and 

• Have a constant voltage output as the rotor speed changes (generator). 

Several types of motor/generator concepts were considered for the spokeless ring flywheel 
system. These are: 

• Annular momentum control device type (reference 17), 

• Ironless armature*. 


*Fisher, R.L., “Ironless Armature Torque Motor,” Contract NAS5-11481, Sperry Marine Systems Division, 
Charlottesville, Virginia, 1968. 
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• Variable strength field (having both stationary armature and field windings), 

• Homopolar d.c. motor. 

None of these were satisfactory by themselves, and the final design is a combination of the 
best features' of the above systems and the results of a number of Goddard Space Flight Center 
programs dating back to the mid 1960s. 

The cross section configuration of the ring rotor motor/generator is shown in figure 6. This 
geometry results in a variable field homopolar d.c. motor (with stationary armature and field 
windings) and will provide a constant voltage variable frequency output from a rotor whose 
speed changes by a factor of 2 (from ^ to co „,„/2). The stationary armature and field 
magnet is continuous around the inside circumference of the ring rotor. The motor/generator 
field is established by two permanent magnet rings (which provide the magnetic bias flux) 
and is then modulated by field coils on the two legs of the stationary magnetic field circuit. 
The stationary armature winding projects into the air gap slot of the rotor moving field. Note 
that the rotor-moving field poles are all north on the upper rings and south on the lower. 
(There are a finite number of equally spaced pole projections on both the upper and lower 
moving iron pole rings.) This configuration (homopolar) means that the flux does not alter- 
nate in direction from pole to pole; rather, salient poles form regions of high flux density, and 
the current return path is in a region of low flux density rather than reversed polarity. 


INSIDE 



Figure 6. Motor/generator system. 
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The field flux included in the motor/generator rotor poles may be utilized as part of the 
magnetic suspension system (see previous section). However, for the sake of simplicity in 
this report, the functions of these systems are considered separately. 

To achieve the proper phase output from the motor /generator the armature windings will be 
electronically switched (commutated) thus providing the most efficient electrical-to- 
mechanical energy conversion and allowing complete symmetry in both motor and generator 
modes. Shown in appendix G are the detailed design equations for the motor/generator 
of the 1 0-kW h rotor described in table 6. The results of these calculations are summarized in 
the list below. 


Homopolar PM d.c. Motor/Generator Characteristics 


Peak Power 
Diameter 
No. Poles 
Flux Density 
Armature Resistance 


= 10 kW 
= 1.8 m (72 in.) 

= 84 (0.81 X 2.28 cm (0.36 
= 0.3 ±0.1 tesla (19.5 ± 6.5 
= 0.087 ohms, R 


160 turns/phase 
30 parallel paths 


X 0.9 in.)) 
kl/in.^ ) 


28 A.W.G. 


Inductance = 3 p henry 

Pole Freq. = 1 1 ,424/s maximum 

• Field Current = 0 

=0.201 V/rad/s 
= 1.07 kg m/A (0.148 lb ft/A) 

= 1 29 V@ 1 02 RPS (75%) 

• Field Current = +3.5 A 

= 0.269 V/rad/s 

K, = 1 .44 kg m/A (0.200 lb ft/A) 

= 115V@68 RPS(50%) 

• Field Current = -3.5 A 

=0.152 V/rad/s 
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Kj = 0.8 1 kg m/A (0. 1 1 2 lb ft/A) 

=130V@136rps 

• Armature Current = 90 A maximum rated 

2r = 705 W ©peak load 

• Constant voltage, charge/discharge 
50% to 100% speed 

1 hr charge to 7 kW h 

IR Drop = 7.8 V at maximum rated current 
Efficiency, at peak load = 93% If = 0 

= 89.5% If = 3.5 

At % CO max speed 
Ij ^ R losses = 420 W 

The electronic commutation function has been successfully applied to a 223-N (50-lbf) 
annular momentum control device (AMCD) ring rotor (reference 17) and the sensors and 
commutating circuitry will be applicable in this design. The output of the motor/ 
generator (see figure 3) will be rectified and inverted to 60 Hz, 1 10 V a.c. power by conven- 
tional means, and no difficulty is expected in this area. 

It is anticipated that the parasitic wattage losses for the motor/generator system will be 
low. The motor/generator is of essentially “ironless” construction meaning that the rotating 
field gap flux does not pass through any stationary magnetic iron; thus, there are no magnetic 
losses due to hysteresis or eddy currents in the armature. Under no-power conditions, the 
only losses are due to eddy currents of the armature wires which reside in the field gap. 
However, these losses may be minimized by reducing the physical size of the wires in the gap 
by means of parallel paths. Since the rotor will change speed over a long period of time 
(typically twice a day), the hysteresis losses of the iron is negligible. The other hysteresis 
losses (the same as those in the magnetic suspension) are due to the mechanical tolerance in 
the shape of the ring and perturbing forces which may cause periodic variations in the moving 
iron ring into and out of the magnetic flux. Since it is doubtful that either of these factors 
will be large, the hysteresis losses are considered negligible. The only other parasitic loss is 
the power dissipation in the field coils. Because the majority of field flux is provided by the 
permanent magnets, this loss can be kept very low (less than 5 percent of the total stored 
energy in the system). 


24 



When power is being drawn from the system, there is an additional loss in the wires of 
the armature. Worst-case calculations have shown that this loss can be easily held to a 
maximum of 5 percent of the total stored energy in the system. 

Efficiency for the charge and discharge cycles is expected to be about 80 percent, and the 
coast (storage) cycle efficiency is expected to be about 90 percent. 

CONCLUSIONS 

There is a demonstrated need for energy storage systems for peak shaving in the electric 
power industry. Energy storage in rotating flywheels provides one solution provided that 
the system life is greater than 20 years, and high energy conversion efficiencies (greater 
than 60 percent) can be achieved. 

One of the most important measures of performance of flywheel energy storage is energy 
density (stored energy per unit of rotor weight). The highest energy densities (and, in a 
general sense, the lowest capital cost of stored energy per dollar) of flywheel rotors are 
obtainable using composite materials in uniaxially stressed geometries. A thin spokeless 
composite ring rotor (ID/OD = 0.9) has a larger energy density than a circular brush superfly- 
wheel and nearly as large (92 percent) as the theoretical maximum for a thin-ring geometry. 

Thus a thin-ring spokeless geometry is extremely favorable for the rotor of an energy storage 
system. 

It has been suggested in this report that attention must be given to all the modules of a com- 
plete energy storage system if long life and high efficiency are to be achieved. In particular 
a general storage system consists of storage, environment, suspension, and input/output modules. 

To best utilize the thin-ring, spokeless geometry as the storage module, a permanent mag- 
netic flux, biased, magnetic levitation system has been configured. This system provides for 
active control in the radial direction with passive magnetic constraints in the axial direction. 
Based on the proposed configuration, little derating of the flywheel energy density is expected, 
and parasitic losses have been shown to be negligible. 

The spokeless ring is also used as the moving element of a specially developed variable field 
homopolar d.c. motor/generator. All parasitic losses are shown to be small, and the system 
functions equally well as either a motor or a generator. 

Design formulas for the parasitic windage loss of the spokeless ring flywheel have been 
developed. Because operation at reduced pressures is a requirement in minimizing windage 
losses, it has been shown that it is better to consider the fluid (air) as a molecular stream 
rather than a simple laminar flowing fluid. 

While this report addresses the application for electric power industry use, it is clear that this 
technology is also compatible with solar- and wind-powered systems. In addition, the storage 
concept is also useful for emergency standby use for hospitals, elevators, computers, and other 
emergency type applications where power interruptions would cause hazards or hardships. 
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SUGGESTIONS FOR FUTURE WORK 


If a prototype of the spokeless, magnetically levitated ring is to produce near-theoretical 
energy densities, some additional theoretical work must be done. In particular, the ring 
element, composed of both composite material and magnetically soft iron, must be optimized 
for mechanical strength, practical fabrication feasibility, and minimum reluctance. The 
following suggestions are offered: 

• The ring design, to be modeled by a NASTRAN finite element program, should 
take into account the anisotropic properties of the composite material and 

the presence of iron in the ring. Various configurations can be investigated 
theoretically and the one with the highest energy density selected. 

• Dynamic analysis of the mechanical/structural system, electrodynamics, and, 
in particular, cross coupling effects should be undertaken. This would include 
coupling of in-plane and out-of-plane rim modes due to radial and axial 
suspension inputs. 

• A cost analysis should be performed on the complete storage system. The best 
way to do this is to take the ring geometry obtained in the first section and 
develop a general computer program which takes into account the total costs 
of all system components. This program should include both projected fabrica- 
tion and material costs and should be capable of determining the best system 
size for a given total energy storage requirement. 

After completion of the analysis, design, build, and test a total mechanical capacitor (7 kW h) 
system. 


Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland November 1975 
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APPENDIX A 


SUMMARY OF UNIT CONVERSIONS 


Given 

Multiolv Bv 

To Obtain 

Specific energy in (N m)/N 

1 

Specific energy in J/N 

Specific energy in in. Ibf/lbf 

3.139 X 10-® 

Specific energy in Wh/lbf 

Specific energy in in. Ibf/lbf 

2.54 X 10-2 

Specific energy in J/N 

Specific energy in in. Ibf/lbf 

6.921 X 10-5 

Specific energy in Wh/kgf 

Specific energy in kWh/lbf 

8.09 X 10® 

Specific energy in J/N 

Specific energy in kWh/lbf 

2205 

Specific energy in Wh/kgf 

Specific energy in Wh/kgf 

367 

Specific energy in J/N 

Specific energy in Wh/lbf 

809 

Specific energy in J/N 

Specific energy in Wh/N 

9.81 

Specific energy in Wh/kgf 

Kinetic energy in in. Ibf 

0.113 

Kinetic energy in J 

Momentum in ft Ibf s 

1.36 

Momentum in J s 

Mass moment in ft Ibf s^ 

1.36 

Mass moment in J s^ 
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APPENDIX B 


EVALUATION OF VARIOUS ROTOR SHAPES 
CONSTANT STRESS DISK (Reference 22) 


Geometrical Configuration 



E 


W 


1 .000 ^ 

T 


a 


max 


E 


w 


W 

V 


a = a . = a constant 
r o 


KsX 


a 

7 


1.000 


1.000 G 

E 

^ = 1.000 7 


2 2 
puj r 


= re- 20 


ID =0 solid 
OD = 


PIERCED DISK (Reference 22) 

Geometrical Configuration 

Both radial and circumferential stresses are considered. 

a = inner radius 

b = outer radius 
T = thickness - constant 




Formulas 


K.E. 


E 


w 


O 

max 


E 


w 


K 


s 


W 

V 

s 


HOOP 


V^poP" rn 

X {b^ + a^} 
7 


(3 + v) 
4 



pco^b^ 


K X 


s 


a 

7 



Geometrical Configuration 

A thin pierced disk where it is assumed = 0 and = po;^b^ (i.e., a/b approaches 1 .0). 



a = inner radius 
b = outer radius 
r = thickness - constant 
~ radial stress = 0 
= tangential stress 


Formulas 

K.E. 

E 

w 

O 

max 


Vapco^ ttt | b ^ 

x|b^ + 

7 

pco^b^ 
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E 


w 


K 


s 



4 


K.x,x(,-(jy) 




NOTE: A given volume could be occupied by a number of independent hoops with 

the effect of the overall a/b being much less than L For this case, the hoop 
formulas apply. (This was the analysis used by Rabenhorst (reference 1 1).) 

SINGLE FILAMENT-ROUND 


Geometrical Configuration 


2b 



Formulas 


2b = length of rod 

R = radius of rod 
Vol = volume of rod 


K.E. 


E 


w 


a (@ r = 0) = 

max'^ ^ 


K 


W 


'/i X pop- X 1/12 |3R^ + 4b^} X Vol 
1/24-^iji X +4b^} 

P0p-\p- 


K X 


s 


a 

7^ 


^ XRXa = K X ax/S-l 
b ^ ( b/ 
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SINGLE FILAMENT-SQUARE 


Geometrical Configuration 



2b = length of rod 

r = thickness of side 
Vol = volume of rod 


Formulas 

K.E. 


a (@v = 0) 

max ^ 


K 


W 


H X X 1/12|t^ + 4b^} X Vol 
1/24£^^|t 2 + 4b^|- 


7 

14 pu3^ b^ 


K X 2. 
s ^ 


1/3 






Ew 


BRUSH-n SINGLE FILAMENTS 


Geometrical Configuration 

It is assumed that n-square or round rods are rotating about a common axis. (See Rabenhorst, 
reference 8 or 13). 
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Formulas 


NOTE: 


K.E. 

E 


a (@ r = 0) 

max ' 


n X (K.E. of 1 square or round rod) 
same as for 1 square or round rod 

same as for 1 square or round rod 


K 


K X2. 

S ^ 


same as for 1 square or round rod 
n X (Ey for 1 square or round rod) 

E / 

— = n X for 1 square or round rod\ 

w ' s / 


Rabenhorst (reference 13) used a total of 20 square rods each 1.27 cm (0.5 in.) 
on a side and 2.25 m (88 in.) total length (2b) for this case: 


n = 20 
K. = 1/3 

and the percent of the swept volume occupied by rods is approximately 15 percent. 
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APPENDIX C 


ID/OD RATIO TO ELIMINATE RADIAL DELAMINATION DUE 
TO CENTRIFUGAL STRESSES 


PROBLEM 

Given a pierced disk where: 

a = inner radius 
b = outer radius 

locale the radius ratio (a/b) where the following condition holds; 


(max) 


< 0.005 


o 


max 


SOLUTION 

Assume the formulas for an isotropic pierced disk apply (reference 22) 


(l+(a/b)2 -l2/b)i . 

^ ^ (ylhf 


{ 3±ii\ (i + . + £a/]^ 

\ 8 ; i ^ ' (j,/b)2 



Where: 


p = mass density 

V - Poisson ratio of contraction in the 
radial direction due to extension in 
the circumferential direction, 0.33 
assumed. 

CO = angular velocity 
= radial stress 
= circumferential stress 


Unit 

kg/m" 


rad/s 

Pa 

Pa 


For convenience, the nondimensional radial and tangential stress are shown plotted versus 
radius ratio (for various a/b ratios) on figure C-1. 

The maximum stress occurs at the inner radius and is given by: 

(^) ^ ^ 

(^)max = poj^b^ (1 - (a/b))^ (C-2) 

Dividing equation C-1 by C-2 and applying the constraint given in the problem statement: 

J2L= (1 - a/b)^ ^ nns 

which can be solved to obtain : 

a /b s 0.9 
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APPENDIX D 


SIMPLIFIED DESIGN EQUATIONS FOR SIZING A STORAGE ROTOR 
PROCEDURE 

The procedure outlined below is applicable to sizing a composite energy storage ring. It is 
assumed that the entire rotating weight is made of Kevlar-49 composite material, and there 
is no derating necessary because of the iron (magnetic flux path elements) which is carried 
on the inside radius of the ring. First 

• Decide on the amount of energy which is to be stored (TE = total stored 
energy) in units of joules (kW h), and 

• Determine the required weight (W) of the ring. 

E =K X 2. 

W S ^ 

where: = shape factor = 0.46 for a thin ring of a/b = 0.9 

S. =0.717.10* J/N (89 Wh/lbf) 

7 

for Kevlar-49 composite being used at 59 percent of its ultimate tensile strength. 

=3.30.10'^ J/N(40.9Wh/Ibf) 

yy _ TE (joules) (N in newtons) 

3.30 X 10"* 

W = (W in IbO 

40.9 

Determine the size of the ring. 

W = ym (b^ - a^) 

= 7T7Tb^ (1 -(a/b)^) 

= 0.597 71-b^ (D-1) 

(Recall a/b is fixed at 0.9 by radial delaminating stress limits.) 
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Choose a reasonable b and solve equation D-1 for t at this point. 

b = known (designer’s choice) 

T = known (from equation D-3 since weight - W - is known) 
a = 0.9b = known 

Determine the speed of rotation by solving equation C-1 for co (assume 0.33 for i^). 


CJ 


■i- 


4 (T, 


dm 


p (3 + v)h^ (1 + [0.81]) 


1 -V 
3 + u 


(D-2) 


u> 




rad/s' 


(D-3) 


If the storage ring is made of Kevlar-49 (see reference 13) operating at 59 percent (this 
represents a safe fatigue limit) of the ultimate tensile strength, then; 

°em ^ 

y = 13.6 k N/m^ (0.05 Ibf/in.^) 


g =9.80 m/s^ (386 in./s^) 

and equation D-3 becomes: 


or 



.. _ 3.37 X 10^^ 

Cu 

b 


rad/s if b in meters 


rad/s if b in inches 


RPM^ 


CO 

2it 


X 60 


rot /min 


Finally, check to see ( 0 ^,)^^^^^^ and (Og)^^^ are all within design specifications. 

EXAMPLE 

Store a total energy of 36 MJ (10 kWh). 

1. TE = 36 MJ (10 kWh) 

2 = 3 6_ X 1 ^ =1091 N ( =244 IbA 

3.3 XIO'^ ' 
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3. 1091 


Assume: 


= 0.597 7rb^ 

7 = 1 3 .6 kN/m^ (0.05 Ibf/in.^) 

b =l,m(39.37in.) 
a = 0 .9 m (35 .43 in.) 

^ = — 1091 = 0.134 m (13.4 cm (5.27 in.)) 

0.597 (136 X 10^)p 


4. 


oj = = M± = 854 rad/s 

b 1 

= IM51X 60 =8163rpm 
2it 


Check the maximum stresses for these conditions using equations (C-1) and (C-2). 
(^r^max a* f = ^b = 0.95 m (37.35 in.) 


(^r)max =(i^)( 13.6X 10 l)(855)^ (1)2 (1 -(0.9))2 
8 9.80 

= 4.22X 10® Pa(612psi) 

(‘2e)max at r = a = 0.9 m (35.43 in.) 


and 


(‘^e^max =( 331) ( 13.6 X lO i ) (855)2 (1^ + 0.2 (0.9)^) 
4 9.80 

= 9.81 X 10® Pa(141,000psi) 


(<^) 

= 4.3 X 10'2 = 0.43 percent 

(‘^.).ax 

which is well within the strength limits of Kevlar-49. 
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APPENDIX E 


WINDAGE LOSSES FOR A RING FLYWHEEL 
WINDAGE LOSSES OF A SPINNING RING 

The following symbols have been used in the deviation which follows: 


Symbol 

Name 

SI Unit 

English Units 

P 

Fluid mass density 

kg/m^ 

Ibf s^ /in.^ 

CO 

Angular velocity 

rad/s 

rad/s 

b 

Outer ring radius 

m 

inches 

a 

Inner ring radius 

m 

inches 

T 

Ring thickness 

m 

inches 

P 

r 

Air pressure 

Pa 

torr 

P 

Power loss 

W 

Whp 

P 

Absolute 

viscosity 

Pa s 

Ibf s/in.2 

Cd 

Drag coefficient 

- 

- 

R 

e 

Reynolds number 

- 

- 


CJT^P 



. 


P 



M 

Mach number 

- 

- 


V tip/V sound 




The drag on a spinning ring (containing no spokes) is skin friction drag which occurs on the in- 
side and outside diametral surfaces and on the top and bottom ring faces. If it is assumed, 
for a worst case condition, that the ring is rotating in stUl air, then the calculation for the 
total drag on the body proceeds as shown. 

NOTE: This procedure must be used (especially at pressures less than 0.13 N/m^ (10"^ torr)) 

since the gas looks more like a molecular stream than a fluid. 

1. Calculate the R^ for both the OD and ID of the ring. Use: 



fx = 2.58 X 10"^ (Ibf s/in}) or 1 .78 X 10"^ (Pa s) for air. This is independent of pressure. 

P 

p=^ 1.15 X 10"'^ (lbfs^/in.^)or(Pj./1.013X 10^) X (1 .23) kg/m^ 

760 

CO = angular velocity (rad/s) 
r = radius (inches) or (m) 

2. Calculate the Mach number for both the OD and ID of the ring. 

V 
TIP 

M= 

V ^ 

sound 

%ound “ room temperature 

3. Calculate the nondimensional parameter for both OD and ID 

Z 

' 1.5M 

4. Use Zj to compute a drag coefficient for both the outside and inside diameters 
by using the following formula (reference 30); 

1-33 p 2Z 

Cn= (EXP(Z2) ERFCZ, -1 + (E-1) 

Z^M L J 

(or the graph shown on page 290 of reference 30) subject to the following 
conditions. 

(a) for Zj <1.0 you are in the free molecular flow region 

c = 0.674 
M 

(b) For 1 < Zj <10 use the graph or formula to obtain 

(c) For Zj > 10 look at the number on the OD. If < 10“^ use the 
laminar flow equations for drag. If the R^ < 1 0^ , use the turbulent 
flow equations. If 10^ < R^ < 10® average the power losses computed 
using both the laminar and turbulent equations. These losses will be 
generally much larger than those computed in (a) and (b). 

5. Compute the power losses according to the formulas shown below. These 
were derived by considering the drag coefficient to; 

• Be a constant for Z^ < 10 

• Be equal to (reference 30, page 78) - 2 > lo and R < in'* 

• Be equal to (reference 30); C„ = Z^ > 10 10® <R^ < 10 
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EQUATION 

Region 

Zj < 10 P (watts) = Cp pco^ 1 .24 (b® - a®) + 3.14 r (b^ + a'*) 

Zj > 10 Rg< 10^ P (watts) = p"-5w2-s (8.82 X 10-^)(b'‘ -a^) + (1.76X lO'^) 

T (b^ + a^) 


Zj>10 Rg>10®P(watts)=p°®w^-®(1.10Xlff2)(b‘*-®-a'^-*) + (2.99X 10'^) 

T(b^-^ +a^-®) 

where the following units must be used: 


p 

fluid mass density 

kg/m^ 

CD 

angular velocity 

rad/s 

b 

outside radius 

m 

a 

inside radius 

m 

Region 

<10 

P (watts) = pco^ 0.14 (b^ - a^) + 0.355 r + a"^) 

Z >10 R 

1 e 

<10^* P (watts) = P®-® 

® (1.20- 10'®)(b‘* -a'^) + (2.39X lO’®) 

Z >10 R 

1 e 

>10^ P (watts) =p°-® 

T (b^ + a^) 

co^-®(2.13X 10-^) (b^-® -a^-®) + (5.76X lO^^) 

where the following units must be used 

f (b^® + a^-*) 

P 

fluid mass density 

Ibf s^ /in.'* 

CO 

angular velocity 

rad/s 

b 

outside radius 

in. 

a 

inside radius 

in. 




APPENDIX F 


FLYWHEEL--SUSPENSION 


To suspend all 1-kg (244-lb) flywheel at its inside diameter magnetically with a passive axial 
(vertical), active radial system: 

Rotor weight = 1 1 1 kg (244 lb) 

Assume design for 2g (gravity) capability 

Load = 222 kg (488 lb) 

Rotor inner circumference = 1,8 m 70.8 in. 

X 3.14 X 3.14 

5.65 m 222 in. 

Force must be produced tangentially at the rotor/stator interface. 


= 384.2 N/m 

5.65 m 

The radial force across this interface is required to be about ten times this force,* 

= 21707 newtons 

Assuming a nominal gap flux density of 0.62 tesla. 


F = 


A 

k 


k = 2.5 X 10-® 
P = teslas 
A = m^ 


A = k £- =k 

A = 0.1422 
(Circumference) = 5 .65 m 

Supplied by motor field: 


21201 = 0.1422 m^ 
0.3844 

= 0.0252 m 

0.0061 m 
0.0191 m 


*See Mechanical Technology Inc., Report 75 TR 21, January 1975. 
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For symmetry, this will be divided into two elements each with two active surfaces. The 
face width of each will be 0.476 cm, further subdivided into teeth of 0.1 59 cm width to 
reduce the sag at Ig to 25 to 50 percent of that dimension. 

The overall axial stiffness should be: 

= 14,244 N/cm 

0.076 cm 

The length of magnet required will be: 

jqj _ 6 X 2 X 8g 
0.000126 

^ 0.62X 2X0.1524 
0.000126 

= 1500 AT 

Using Sm COj magnet at H = 3200 AT/cm: 

8 = 0.469 cm 

m 

At this operating point, the magnet can supply 0.4 tesla. Since the required flux across the 
gap is 0.62 tesla at a total width of 2.52 cm and allowing for leakage: 

w =mXM78XM = i.33cm 
■" 0.4 

Since the force is proportional to the square of the flux density, a two-to-one change in 
force can be achieved by a ±25 percent modulation of the magnet ampere turns. The effec- 
tive air gap, for the control winding, includes the permanent magnet; therefore, the coil 
ampere turns required: 

%(l + X ^=201) 1500 AT =462 AT 

V 1.33 0.468 / 

To modulize the assembly, a group of four pole spans will be considered. The core cross 
section will be 36.12 cm^ (5.6 in.^) if wound around the magnet. Assuming a path length 
of 30.5 cm (12 in.) per conductor and maximum of 40 W per coil, 360 turns of No. 27 A.W.G. 
wire would result in a peak coil current of 1 .3 A and a coil resistance of 18.5 ohms. A 
0.65-cm^ (0.1-in.^) coil cross section would be required. 
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APPENDIX G 


MOTOR/GENERATOR CALCULATIONS 
Schematics are shown in figures G-1 through G-4. 

Type: Electronically Commutated Ironless Armature d.c. Variable Field, Induced Flux 
Homopolar. 


Power 

Input/Output 

10,000 W 

Diameter 

Nominal 

1.8 m 

Speed 


641 ±214 rad/s 

Voltage 

Nominal 

120 V 

Current 

Max 

90 A 

No. of Poles 


84 

No. of Phases 


3 

No. of Conductors 
Per Phase 


136 

No. of Parallel 
Paths 


30 

Wire Size 


28 A.W.G. 

Resistance 

delta 

0.087 ohms 


per phase 

0.1 3 ohms 


per path 

3.9 ohms 

Flux Density 

Pole Gap 

0.3 tesla (19.5 kl/in.^) 


Interpole 

0.15 tesla (9.75 kl/in.^) 

Pole Area 


0.0002 m^ (0.324 in.^) 

Pole Span 


0.068 m (2.7 in.) 
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E = N-^ X 10'*V 

dt 


N = 160 conductors 

A0/cm = 1 .38 kilolines/cm (A0/in. = 3.5 kl/in.) 


CO = 102 rps (75% speed) 

E = 129 V 

a 


Force = = Newtons 

= 

8.8 X 10'® lb 

Aj3 = 0.15 tesla 


9.75 kl/in.2 

I = 2.9 A 


2.9 A 

C = 0.0091 m 


0.36 in. 

= 0.0036 N/conduct or 

= 

0.00895 Ib/conductor 

X 160 conductor 


X 160 conductor 

0.63336 N 


0.1433 lb 

X 30 oaths 


X 30 oaths 

19.0 N 


4.3 lb 

X 0.914 m 


X 3 ft 

17.37 Nm 


12.9 lb ft 

X 102 rps 


X 102 ros 


11200W 11217W 


I™,v=87A I2r=659W 

max 

Inductance, from previous design: 

L = 0.00041 X ( 1 2 = 3 X lO'* henry 
\1872/ 

L/R = 3 = 34.4 /is 

0.087 

Frequency, co ^ poles = 11 ,424 Hz 

Pole face size = 0.92 X 2.29 cm (0.36 X 90 in.) 

Conductor turn length = 2 (2.29 cm + 3.43 cm) (2(0.9 in. + 1 .35 in.)) 
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1 1 .44 cm/conductor 
X 2 conductor/pole 

22.88 

X 80 poles (active) 
1830 cm 


4.5 in./conductor 
X 2 conduct or/pole 
9.0 

X 80 poles (active) 
12 I 720 in. 

60 ft 


X 0.00214 ohms/cm (#28 A.’ 

30 I 3.92 ohms/path 
0.1309 ohms/phase 
X 0.67 

0.087 ohms (Delta) 
Total conductor length = 

1646 m = 

X 0.00745 ke/m 
12.3 kg 

Motor-only rotor iron = 645 cm^ 

5.27 kg 


.G.) X 0.0653 ohms/ft (#28 A.W.G.) 

30 I 3.918 ohms/path 
0.1306 ohms/ phase 
X 0.67 

0.087 ohms (Delta) 


X 3 X 30 = 5400 ft 

X 0.0005 Ib/ft 
2.7 lb 

39.4 in.3 
1 1 .6 lb 


FIELD-PERMANENT MAGNET 


Flux per pole = 6.3 kilolines 
Interpole leakage =6.3 kilolines 
Pole gap = 0.635 cm 

Stator gap = 0. 1 524 cm 


@0.2 tesla (@ 19.5 kl/in.^) 
@0.15 tesla (@ 9.75 kl/in.^) 
(0.25 in.) 

(0.06 in.) 


HF =BS— 
3.192 


ampere turns 


= 2452 AT 


- 0.95 cm (0.375 in.) Samarium Cobalt 


FIELD-CONTROL COIL 

Control flux density = 0.1 tesla (± 6.5 kl/in.^) 


HF 

2212 =973 AT 


3 

9 

16.64 kl/4 poles 


1.89 tesla (122 kl/in.^) 

Area = 

0.88 cm^ (0.136 in.2) 
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diameter 

21 coils, 

974 AT 
0.17 A 

Wire size, 36 A.W.G. 


1.05 cm (0.416 in.) 
0.17 A/coO 

5800 turns 


Coil resistance, 700 ohms 

Max dissipation 20 W/coil 

21 

420 watts, field 


FILAMENTARY MOTOR 



Figure G-1 . Rotor/stator cross section. 



ROTOR 

POLES 


ROTOR 



Figure G '2. Rotor/stator top view. 



Figure G-3. Stator/motor suspension 






2.54 cm 
(1.00 in.) 


Figure G4. Rotor. 
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